Subspace Confinement: How Good is your Qubit? 
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The basic operating element of standard quantum computation is the qubit, an isolated two-level system that 
can be accurately controlled, initialized and measured. However, the majority of proposed physical architectures 
for quantum computation are built from systems that contain much more complicated Hilbert space structures. 
Hence, defining a qubit requires the identification of an appropriate controllable two-dimensional sub-system. 
This prompts the obvious question of how well a qubit, thus defined, is confined to this subspace, and whether 
we can experimentally quantify the potential leakage into to states outside the qubit subspace. In this paper we 
demonstrate that subspace leakage can be quantitatively characterized using minimal theoretical assumptions 
by examining the Fourier spectrum of the oscillation experiment. 

PACS numbers: 03.67.Lx, 03.65.Wj 



I. INTRODUCTION 

The issue of subspace confinement for qubit systems is fun- 
damental to the primary operating assumptions of quantum 
processors. The concepts of universality, quantum gate op- 
erations, algorithms, error correction and fault-tolerant com- 
putation hinge on the precept that the fundamental quantum 
system is an isolated, controllable, two-dimensional system 
(qubit). 

It is well known that most of the physical realizations 
of qubits are in fact multi-level quantum systems, which 
can theoretically be confined to a two-dimensional (qubit) 
subspace. Important examples range from super-conducting 
qubits JH 13, 13I] to atomic systems such as cavity-coupled 
color centers 14, IH 0] and ion traps JtII- In the former sys- 
tems, a qubit is generally defined as the subspace (of the 
full Hilbert space) spanned by the two lowest energy states 
in an arbitrarily shaped potential such as the washboard po- 
tential of cuiTent-biased Josephson Junctions |^. How- 
ever, the potential number of valid quantum states within each 
well is not limited to two, and quantum gates, especially if 
sub-optimally implemented, may inadvertently populate other 
confined states. Similarly in ion trap systems, a qubit is usu- 
ally defined by two electronic states of an ion, either two hy- 
perfine levels or a ground state and a meta-stable excited state, 
but once again there exist many other electronic states. Hence 
a more stringent definition of a qubit would consist of a two- 
level quantum system with classical control confined to the 
unitary group SU{2). 

The ability to initialize, operate and measure completely 
within the two-level subspace representing "the qubit" is vi- 
tal to the successful operation of any large scale device con- 
structed from such quantum s^tems. Standard quantum error 
correction protocols (QEC) lfl0[[TTl[l2ll generally assume that 
the qubit system is precisely confined to the two-level sub- 
space and that all quantum gates operate only on the qubit de- 
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grees of freedom. If poor control or environmental influences 
inadvertently results in non-zero population of higher levels, 
leakage correction protocols are necessary. 

The issue of subspace leakage in quantum processing has 
been addressed in depth from the standpoint of error correc- 
tion. Work by Lidar ifisi [Till examined the construction of 
Leakage Reduction Units (LRU's), which use modified puls- 
ing techniques to ensure that any unitary dynamics outside 
the qubit subspace can be compensated for which has been 
adapted specifically for super-conducting svstems |15il. An- 
other type of LRU's uses quantum teleportation lUql to map 
a multi-level quantum state back to a freshly initialized two- 
level qubit. Finally, active detection such as non-demolition 
measurements (which detect population in non-qubit states 
without discriminating between t he q ubit states) can be per- 
formed on the system UtI [isl [l9l l20ll . If an out-of-subspace 
detection event occurs the leaked qubit is re-initialized or re- 
placed. The inclusion of LRU's based on teleportation has 
been investigated within the context of fault-tolerant quantum 
computation f23\ and shows that, in principle, the inclusion of 
leakage protection does not adversely affect large scale con- 
catenated error correction. 

Although these schemes are viable methods to detect and 
correct for improperly confined qubit dynamics, they can be 
cumbersome to implement and many systems admit, in princi- 
ple, sufficiently confined Hamiltonian dynamics so that leak- 
age could be expected to be heavily suppressed. For example, 
for ion-trap qubits controlled by lasers, leakage to other ionic 
states can be made negligible by employing very finely tuned 
lasers and sufficiently long (and possibly optimally tailored) 
control pulses. Advances in qubit engineering may therefore 
allow us to eliminate or at least substantially reduce the need 
for laborious leakage detection/prevention schemes in many 
cases, provided that we can experimentally ascertain suffi- 
ciently high confinement of manufactured qubits under classi- 
cally controlled Hamiltonian dynamics. 

In this paper we present a simple generic protocol to esti- 
mate qubit confinement, or more precisely, establish bounds 
on the subspace leakage rates, for "quality control" purposes. 
The main goal is to allow us to empirically detect inferior 
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qubits by using readily obtainable experimental data to derive 
tight bounds on the subspace leakage of the system. This pro- 
tocol would represent one of the first steps towards full system 
characterization ll22ll2l |24i B . 

Section briefly outUnes the basic assumptions with re- 
spect to the measurement and control model and the motiva- 
tion for the proposed protocol. Section|lIl]discusses the basic 
mathematical properties of qubit oscillation data and shows 
how a minimal amount of information obtained from the oscil- 
lation spectrum can be used to derive empirical bounds on the 
subspace leakage rate, and that these bounds are very tight for 
the high quality qubits required for practical quantum compu- 
tation. In section HVl the effects of finite sampling are consid- 
ered and studied using numerical simulations. SectionlVjcom- 
pares the efficiency of bounding confinement using the pro- 
posed scheme versus alternative approaches such as detection 
of imperfect confinement by identifying additional transition 
peaks within the Rabi spectrum. Finally, section |Vl] briefly 
examines the effects of decoherence. 



II. MOTIVATION AND PRELIMINARIES 

Estimation of qubit confinement represents one of the first 
major steps in full qubit characterization. Therefore, the pro- 
tocol should not be predicated on the availability of sophis- 
ticated measurements or control, and should be amenable to 
automation so that it could be used in conjunction with a po- 
tentially automated qubit manufacturing process. The bounds 
on the subspace leakage will be based on the observable qubit 
evolution under an externally controlled driving Hamiltonian. 
We assume that our classical control switches on the single 
qubit dynamics and that the governing Hamiltonian is piece- 
wise constant in time. Hence the Hamiltonian induces the uni- 
tary operator U = e^'^*, with h ^ 1. 

Although this assumption may not be applicable to all sys- 
tems, e.g., systems subject to ultra-fast tailored control pulses, 
it is not as restrictive as it might appear. It is generally be valid 
for systems such as quantum dots or Josephson junctions sub- 
ject to external potentials created by voltage gates if the gate 
voltages are (approximately) piecewise constant. It is also 
a good approximation for systems subject to time-dependent 
fields such as laser pulses in a regime where the rotating wave 
approximation (RWA) is valid and the pulse envelopes can 
be approximated by square-waves. In this case, the Hamil- 
tonian relevant for our purposes is the (piecewise constant) 
RWA Hamiltonian determined by the amplitudes, detunings 
and possibly phases of the control pulses. This model can even 
be valid for other pulse shapes if the Hamiltonian is taken to 
be an average Hamiltonian describing the effective dynamics 
on a certain time scale (beyond which we do not resolve the 
time-dependent dynamics). However, the main focus of the 
paper is not when the dynamics of a system can be modeled 
in this way, but rather how to assess subspace confinement for 
systems where this model of the dynamics is valid. 

Assuming the effective control-dependent Hamiltonian 
H = H[f] is constant for < t < tk, where / is the clas- 
sical "control knob" parameter, the evolution during this time 



period is given by the unitary operator U{t) = e~*^*. Al- 
though H will generally depend on control inputs, we shall 
omit this dependence in the following for notational conve- 
nience. The driven system generally undergoes coherent os- 
cillations, which are often referred to as Rabi oscillations, es- 
pecially for optically driven systems in the RWA regime. Al- 
though our model is not limited to these systems, we shall 
use the terms coherent oscillations and Rabi oscillations inter- 
changeably throughout this paper. 

The measurement model assumed is crucial to the relevance 
of the protocol. Some standard measurement models in quan- 
tum computation assume the ability to detect both the |0) and 
|1) states independently (such as SET detectors in solid state 
designs ll26ll27il28ll ). In this case, estimating subspace leak- 
age is fairly straightforward and requires only repeated mea- 
surement of the system while undergoing evolution. The leak- 
age is simply given by the deviation of the cumulative proba- 
bility of measuring |1) or |0) from unity. However, this mea- 
surement model is not realistic for the majority of proposed 
systems. 

Color centers and ionic qubits use externally pumped tran- 
sitions to discriminate between a light state (= |0)) and 
other "dark" states, while readout in super-conducting sys- 
tems ll29i [30I1 involves lowering a potential barrier such that 
only one of the qubit states can leak to an external detection 
circuit. The measurement outcome of the indirectly probed 
state is inferred from the non-detection of the directly mea- 
sured state and for such measurement models estimating con- 
finement is more complicated. Hence this paper utilizes the 
latter model in order to quantify confinement. It should be 
noted that we are not considering the concept of weak mea- 
surement, in each case we assume that the measurement of 
the system causes a full POVM collapse of the wavefunction. 
We also assume that the measurement apparatus has been suf- 
ficiently characterized. In order to to successfully implement 
computation, readout fidelity should ideally be of the same or- 
der as general systematic and decoherence errors. Therefore, 
characterization is initially required to ascertain the error rate 
associated with measurement which can then be incorporated 
into calculations of confinement. 

Strong non-qubit transitions can still be identified directly 
via modulations in the Rabi oscillations data as shown in 
Fig-ID' for a three-state system evolving under the trial Hamil- 
tonian 

/ 1 0.5\ 
H,n=\ 1 1 . (1) 
\0.5 1.5 J 

However, the Rabi oscillation data for the modified three-state 
Hamiltonian, 

/ 1 0.0l\ 
Hn= \ 1 1 \ , (2) 
\0.01 1.5 J 

depicted in Fig.[T^ shows that an apparent lack of modulations 
in the Rabi oscillation data is not proof of perfect confine- 
ment, and that quantitative measures of confinement or sub- 
space leakage and experimental protocols are needed. 
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m. ESTIMATION OF SUBSPACE LEAKAGE 

By defining the projection operator onto a two dimensional 
subspace, 11 ~ |0)(0| + |1)(1|, subspace leakage is given by, 

e=l-Tr[np], (3) 



with p = U'''{t)\0){0\U{t). Unfortunately, we cannot calcu- 
late e directly without knowledge of the Hamiltonian. How- 
ever, we can estimate subspace leakage experimentally from 
standard Rabi oscillation data. 



a) b) 




50 



time (Inverse Energy) time (Inverse Energy) 



FIG. 1: Modulations in the Rabi oscillations of a three-level system driven by the Hamiltonians, Hm and Hn- Fig. b) provides clear evidence 
that this system is not a qubit, while Fig. a) appears to show perfect confinement. However the analysis in the following sections will show 
that the subspace confinement for the system in Fig. a) is also not sufficient for large-scale QIP applications. 



A. Perfect confinement 

Consider a general A'^-level system undergoing coherent 
evolution via a driving Hamiltonian in the closed system 
case of no environmental decoherence. If confinement under 
this Hamiltonian is perfect, Hn has a direct sum decomposi- 
tion, 

Hn = H2x2® H(^n-2)x{N-2) (4) 

where iJ2x2 represents the control Hamiltonian confined to 
the qubit subspace, span{|0), |1)}, the state |0) being defined 
by the measurement, and the excited state |1) by the allowed 
transition. For our measurement model the observed Rabi os- 
cillations have the functional form f{t) = \{0\UN{t)\0)\'^ ■ 
As there is no coupling between the state |0) and states out- 
side the H2x2 subspace, we can expand f{t) by diagonalizing 

U2x2{t) = CXp(-ii?2x20 

fit) = |(0|^tdiag{e-'^°*,e-*^i*}A|0)|2 

= ||co|'e-*^°* + |ci|2e-»^^*|2 (5) 

= |co|" + |ci|" + |concip(e^-oi* + e-'"-*) 

where U2x2{t) = A'<e-'"'''A, A\0) = c„\0) + ci|l), ujqi = 
Ao — Ai and {Aj} are the eigenvalues of H2x2- For perfect 



confinement, H2x2 induces coherent oscillations between the 
two qubit levels at a Rabi frequency given by the difference in 
the eigenvalues. Taking the Fourier transform of f{t) gives 

FH=FT[/(t)] = (|co|^ + |cin<5M 

+ \co\^\ci\^6{uj^ujoi) + \co\^\ci\^5{uj + ojai). 

Conservation of probability (total population) thus implies 

(IcoP + lciH^ ^ |co|4 + |ci|4+2|conci|2 = l,andhencethe 
heights of the two Fourier peaks for perfect confinement will 
satisfy the relation ho + 2/io,i = 1, where ho = |co|^ + |ci|^ 
and /io,i = |cop|cip. 



B. Imperfect confinement 

If the system experiences leakage to states outside the qubit 
subspace then the corresponding control Hamiltonian Hn can 
no longer be reduced to a direct sum representation (|4]i but 
it can be diagonalized H^ = diag[{Aj}], {Aj} being the 
eigenvalues of Hn, and the propagator UN{t) expressed as 
UN{t) " A''e^^^'^*A. The Rabi data is now a linear superpo- 
sition of multiple oscillations corresponding to different tran- 
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sitions of the A'^-level system 

fit) = \ {0\A^e-'"^*A\0)\^ 



N-l 

E 

a=0 



(7) 



a.b 



and the corresponding peak heights in the Fourier spectrum 
can be expressed in terms of the expansion co-efficients, 

= J2a=oCa\a), as, 



ho=^\Ca\'^, ha,b = \ca\'^\cb\'^, a ^ b. 
Conservation of probabihty leads to 

.N-l .2 N-l 

a#6 



(8) 



2| |2 



(9) 



Imperfect confinement implies /lo + 2/io,i < 1. We see from 
this analysis that the subspace leakage e is determined by 
the cumulative amplitudes of all non-qubit states for a given 
eigenstate of i^^r, which can be calculated from all the peak 
heights in the Fourier spectrum. 



E 



hh.c 



b,c^ a. 



(10) 



However, exact calculation of e requires identification of all 
peaks in the Fourier spectrum and knowledge of which peak 
corresponds to each trasition. It is therefore desirable to de- 
rive bounds on the subspace leakage that only involve a few 
dominant and thus easily identifiable Fourier peaks. 



C. Bounds on subspace leakage 

We can derive upper and lower bounds on e using only the 
heights of the primary spectral peaks ha and /io,i. 



a^n.l 

^il-ef+Y \Ca\\ 



Provided J^a^a i 1' subspace leakage is reason- 

ably small, we obtain a tight lower bound for e as a function 
of only the two major peak heights: 



.-. e > 1 - y/ho + 2ho,i. 



(12) 



The upper bound for e can also be calculated quite easily. 
Recall that 



E 



E - E i^^i 
> E i^^i'- 

Comparison with ( fTTT i thus immediately yields 



to 



2/io,i < (1 - e)2 + £2 ^ 1 - 2e + 2e^ 



which can be solved for e 
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£<-(!- v/2/io + 4Vi - 1) 



(13) 
(14) 

(15) 



The other solution to Eq. (fl4] i is invalid as a bound due to the 
asymptotic behavior of both the upper and lower bound 



lim min(e) = 0, 

{ho+2ho,i)^l 

lim max(e) = 0. 

(ho+2/io,l)^l 



(16) 



Since the second term in ( fT3b represents the heights of all the 
Fourier peaks «of associated with the |0) ^ |0) ^ \a) 
or |1) ^ \a) transitions, for \a) ^ For a well confined 
system this is a very small correction to e^, consequently the 
bound is again strong. 

Therefore, the subspace leakage e is bounded above and 
below by 

1 



1 - ^ho + 2ho,i <€<-{!- ^2ho + 4/io,i - 1). (17) 

Note that this double inequality involves only the two main 
peaks in the Fourier spectrum, i.e., we can bound the subspace 
leakage without determining the heights of all peaks. 

For the trial Hamiltonians ([T]i and Q we obtain the follow- 
ing bounds 



0.0497 <eH„ < 0.0511, 
3.9754 X 10"* <eH„ < 3.9762 x 10"^ 



(18) 



while the actual values of and e^^ are 

e/f„ = 5.11 X 10-2, e^,. = 3.9762 x 10"^ (19) 

In both cases the upper bound for e equals the actual value of 
e. This is due to the fact that both systems are of dimension 
three, and when estimating max(e) we neglected terms of the 
form 



E 

{a,b)^(0,l),ayib 



Ca\W, 



(20) 



which naturally vanish for a three-level system. 

Fig. [2] shows how the bounds ( [TT] ) for e converge as con- 
finement increases (7 — > 0) for the test Hamiltonian, 



Hi 



^0177 
110 
7 1.5 

V7 l.7y 



(21) 
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0.14 



FIG. 2: Upper and lower bounds on e for the four-level trial system 
governed by the Hamiltonian l l21b . characterized by a static coupling 
between the qubit states and a variable coupling 7 to two higher lev- 
els. As 7 — > the subspace leakage approaches and the bounds for 
e become more accurate. 



IV. FINITE SAMPLING FOURIER ANALYSIS 

The previous section details how quantitative bounds on 
the subspace leakage can be obtained, in principle, from the 
Fourier spectrum of the Rabi data. However, to translate this 
method into a viable experimental protocol we need to con- 
sider the effects of finite sampling and taking the discrete 
Fourier transform (DFT), which raises several issues. 

First the Nyquist criterion for sampling ifsill must be satis- 
fied, i.e., to avoid aliasing, some rough estimate of the Rabi 
period TRabi is needed to guarantee that at least two sample 
points are chosen per oscillation period, i.e., Ai < TRabi/2. 
The second issue that must be considered is the resolution of 
the Fourier spectrum. The frequency resolution Aw is given 
by Acj = 27r/<oh, with tob the total observation time of the 
Rabi signal. If the control Hamiltonian induces a non-qubit 
transition with a frequency within Alu of the primary peak 
then the DFT will combine the amplitudes for qubit and non- 
qubit transitions in the same frequency channel thus leading 
to an overestimate of ho.i and hence qubit confinement. To 
avoid such problems it is necessary to ensure that the total ob- 
servation time tob is long enough. Thus, some estimates of the 
system parameters are required, although these do not need to 
be very accurate and will generally be known on theoretical 
or experimental grounds. 

Finally, the DFT has the property that a pure sinusoidal sig- 
nal will approach a delta function if there is zero phase differ- 
ence between the start and the end of the observed signal. If 
this phase matching condition is not met then all frequency 
peaks will broaden. Phase matching for system identifica- 
tion has already been addressed for the identification of sin- 
gle qubit control Hamiltonians in Izsll and we will follow the 



same approach, which essentially involves truncating the Rabi 
oscillation data at progressively greater values of tob such as 
to maximize the trial function 



P(U) 



2F{up) - Fjup - 1) - Fjup + 1) 

FiLUp-l)+F{LOp + l) 



(22) 



where F{uj) represents the amplitude of the Fourier Spectrum 
at frequency uj and ujp represents the frequency of the max- 
imum Fourier peak. The value of tob where P{tob) is maxi- 
mized represents the cut off time to the Rabi signal that pro- 
duced the best phase matching for the DFT. 

To simulate real experiments we numerically propagate 
the initial state |0), under the Hamiltonian H, by U{tk) = 
exp(— lifei?) for discrete times tk ~ kAt where k = 
0,1,..., if and KAt = tob- A single measurement at 
time tk is simulated by mapping the target state U{tk)\Q) to 
{|0), |1)}, where the probability of obtaining is given by 
Pq — |(0|?7(tfc)|0)p; the ensemble average at a single time 
tk is determined by dividing the number of zero results by 
the total number of repeat experiments A^e • For the following 
numerical simulations we shall use the trial Hamiltonians 



/O 1 
'110 



and 



Hb = 



0\ 

' 

1.5 

1.7 

Vo 02/ 



/ 1 0.01 0.005 0\ 

110 00' 

0.01 1.5 

0.005 1.7 

V 02/ 



(23) 



(24) 



where Ha represents a five-level system with a perfectly de- 
coupled two-level subspace consisting of the two lowest en- 
ergy states, while Hb represents a five-level system with weak 
coupling between the qubit sub-manifold and two of the up- 
per levels. We only consider Hamiltonians that have couplings 
between the |0) state and higher levels, as this state is fixed by 
the measurement basis. We are therefore free to diagonalize 
the lower block of the Hamiltonian, which also helps to sim- 
plify the comparison between different systems. 

The out-of-subspace coupling in Hb was chosen such that 
the leakage from the qubit subspace e « 7 x 10^^ is small (too 
small to cause noticeable modulations in the Rabi oscillations) 
yet significant (in fact above certain critical thresholds) for 
quantum computing applications. The part of the Hamiltonian 
governing the qubit dynamics was chosen arbitrarily and is 
common to all the Hamiltonians examined within this paper 
to maintain consistency between different simulations. The 
accuracy of the protocol is not affected by the choice of single 
qubit dynamics. 



A. Estimating uncertainty in leakage bounds 

Estimating uncertainties in the bounds for e is crucial since 
for the majority of qubit systems it will be practically impos- 
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sible to prove that the evolution of the system under a given 
Hamiltonian is completely confined to the SU (2) subspace, 
i.e., e = 0. Instead, in practice it is sufficient for quality con- 
trol purposes to experimentally confirm that the leakage from 
the qubit subspace is below a threshold value where it can ef- 
fectively be ignored, i.e., it is the upper bound max(e) that is 
relevant. The accuracy of our estimate for max(e) will be pri- 
marily limited by our ability to accurately determine the main 
peak heights and Hqj due to projection noise induced by 
the DFT. 

Quantifying this uncertainty is relatively straightforward. 
Defining the noise function z^(cj) of the Fourier spectrum 
to be the amplitude i'(ti') of each Fourier channel excluding 
ho = ^(0) and ho.i = F{u)p), the uncertainty in /lo and 
/io,i is given by the standard deviation of the noise function 
5h = sd[v{Lo)]. From this we can derive the uncertainty asso- 
ciated with max(e) = e„. 



dho 



\dho 



(ShoY 
( <9e„ 



01 



35h 



]Sho6hri 1 = — , 
01 J ' 2V2/10 + 4/11-1 

(25) 

Seu can be reduced by increasing the number of ensemble 
measurements taken at each point in the Rabi cycle. Fig- 
ures |3] and 2] show how the estimate for converges as A^e 
is increased for the Hamiltonians ( |23] ) and ( l24l i. respectively. 
It should be noted that e„ > 0, hence for each plot the lower 
error bars should only extend to the zero point, but keeping 
the error bars symmetrical around the data point makes the 
convergence behavior clearer. For large values of Ne, e„ con- 
verges to zero for the perfectly confined system governed by 
Ha but the non-zero value w 7 x 10~* for the imperfectly 
confined system described by Hi,. The respective observation 
times for each Hamiltonian were chosen to be tot ~ SOTRabi 
to ensure that all peaks are resolved, i.e., there are no contri- 
butions from additional transitions present within Aoj of the 
primary peak. 



B. Numerical tests of error bound accuracy 

To test the overall accuracy of the uncertainty estimates 
for e„ we can expect to obtain from realistic Rabi oscillation 
data, we calculated the distance between the simulated value, 
eu, and the analytical value, e^, calculated directly from the 
Hamiltonian using Eq.fTSlas. 



diH,) ^ \e.aiH,) - e'jHk 



(26) 



where k E [a,b] and 5d{Hk) is the error in d resulting 
from the error associated with estimating e„(i/fe). We first 
calculated the distance d{Hk) and Sd{Hk) for 5000 simu- 
lated runs of two known trial Hamiltonians (Ha and Hi,) with 
e'uiHa) = and e'^{Hi,) ~ 7 x IQ-"*, respectively. The dis- 
tributions of d{Hk) for Ha and Hb (with Ne = 1024 and 
tob as in Figs [3] and |4] are shown in Figs [5] and |6] respec- 
tively. The average error 3Sd{Hi), i e [a,b], was given by 



FIG. 3: Convergence of eu as the number of ensemble measurements, 
TVe, is increased for a system governed by the Hamiltonian J24t . char- 
acterized by perfect subspace confinement. The solid line represents 
the actual value of eu{Ha) ~ 0. Note error bars should only extend 
to zero as tu{Ha) > 0. 



— 1- 



FIG. 4: Convergence of as the number of ensemble measurements, 
iVe, is increased for the imperfectly confined system governed by 
the Hamiltonian i23t . The solid line represents the actual value of 



eu{Hb) ~ 7 X 10- 
euiHi) > 0. 



Note error bars should only extend to zero as 



i5d{Ha) 



4.92 X 10"'', encompassing 99.9% of the data, 
I w 5.02 X 10-4, encompassing 99.8% of the 



and 36d{Hb) 
data, respectively. 

Next we examined how the protocol behaves when simulat- 
ing a large number of randomly selected multi-level Hamilto- 
nians. For these simulations we choose A^-level Hamiltonians 
of the form 

9 9 

Hn = ^Sfc|A:)(fc| + |0)(1| +^afc|0)(A;| +h.c. (27) 



k=0 



k=2 



with {Ek} = {0, 1, 1.5, 2, 2.4, 2.5, 2.9, 3, 3.3, 4}. The vector 
a = [a2, . . . , flg] was then chosen at random in two stages. 
First the dimensionality of a is randomly selected, allowing 
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d{Ha) 



FIG. 5: Distribution of d{Ha) for 5000 separate simulations. The 
average of the error, 3Sd{Ha) is also shown, with approximately 
99.9% found within Sa of d = 0. 




FIG. 6: Distribution of d{Hb) for 5000 separate simulations. The av- 
erage of the error, 3Sd{Hb) is also shown, with approximately 99.8% 
found within 3a of d = 0. 



the Hamiltonian to coherently drive any multi-level system, 
N £ [2, 3, .., 10]. The non-zero coupling values were then 
randomly assigned such that each element of a was approxi- 
mately two orders of magnitude less than the qubit coupling 
term to ensure that all of the multi-level systems had high con- 
finement. 

We randomly generated 5000 of these Hamiltonians and 
d{Hk) = \tu{Hk) — e'u(-fffe)| was calculated. The average 
(analytical) value of €[XHk) for these 5000 trial Hamiltonians 
was found to be e'^{Hk) = 1.68 x 10"''. We then examined 
the ratio. 



R = 



Num{(d(gfc) - 36d{Hk) < 0)} 
5000 



(28) 



indicating the percentage of successful estimates of the sub- 
space leakage within 3a. This ratio was calculated to be 
R = 99.9%, with the confinement estimates being outside the 



error bounds for only three of the randomly generated Hamil- 
tonians. 

These results are consistent with the expectation that ap- 
proximately 99.7% of the data should lie within 3(t of the 
mean and demonstrates that our methodology for characteriz- 
ing subspace leakage can indeed be expected to yield accurate 
upper bounds on the subspace leakage in the vast majority of 
cases. 



V. EFFICIENCY OF THE PROTOCOL 

The protocol presented in the previous section allows us 
to determine quantitative bounds on the subspace leakage for 
imperfect qubits by determining only the main peaks in the 
Fourier spectrum. An alternative strategy is to try to iden- 
tify all peaks in the Fourier spectrum. The presence of any 
peaks in addition to the two main peaks is indicative of sub- 
space leakage and a quantitative estimate of the leakage rate 
can be obtained by determining the heights of the additional 
peaks. Both approaches have potential advantages and disad- 
vantages. The former approach requires only the identification 
of the two main peaks but these need to be clearly resolved 
and the peak heights determined with high precision. The lat- 
ter approach does not require precise estimates of peak heights 
but relies on the detection of additional peaks, which for high 
confinement will be much smaller than the major peaks, and 
are likely to be difficult to discriminate from the noise floor 
This raises the question which strategy is more efficient to de- 
cide if the subspace leakage for a given qubit is below a certain 
error threshold. 

To answer this question, we performed a series of numeri- 
cal simulations comparing the total number of measurements 
required to ascertain that the lower bound on the leakage rate 
ci = 1 — y^fiQ + 2/10.1 > within error bounds, versus iden- 
tifying a statistically significant third peak in the Fourier spec- 
trum, indicating an out-of-subspace transition, for various trial 
Hamiltonians. For the purpose of the simulations we consider 
the following trial Hamiltonians 




(29) 



representing a system with a variable coupling 7 to a third 
level, as well as the four-level system governed by the Hamil- 
tonian (I2TI 1 and a six-level system governed by 



Ha 



/O 1 7 7 7 7 \ 

1 1 

7 1.5 

7 1.7 

7 1.9 

V7 2.2/ 



(30) 



representing systems with variable but equal coupling to be- 
tween one and four out-of-subspace levels, respectively. 

The lower bound, ej, is taken to be non-zero for a discrete 
data set, if the analytical value e[ of the lower bound calcu- 
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lated directly from the Hamiltonian exceeds six times the un- 
certainty, 5{ei), for the discrete data calculated from the sim- 
ulated Fourier spectrum, i.e.. 



e[ ~ 6S{ei) > 0, 



(31) 



2y^ho + 2ho,l 



Six times the uncertainty in e; represents the total distance 
between the maximum and minimum possible value of e; (us- 
ing a 3cr upper and lower confidence bound) and this interval 
should be smaller than the analytical value, ej . 

A peak F{uj') in the discrete Fourier spectrum is taken to 
be significant if it is more than three standard deviations Sh = 
sd[i'{uj)] above the projection noise floor i.e.. 



F{lu') - D{lu) - 35h > 0. 



(32) 



This definition will underestimate the number of ensemble 
measurements required slightly as it only represents the point 
where the third peak is greater than at least 99.7% of the noise 
channels. 



- - Usmg coiifmement equations 

— Ideiiufying a third peak. 



.002 .0046 .0072 .0102 .0134 

Confinetnent (e) 

FIG. 7: Number of ensemble measurements required to ascertain 
statistically significant subspace leakage (imperfect confinement) for 
the three-level system governed by i29l as a function of the (analyt- 
ically calculated) confinement using the confinement equations l |31| l 
and by directly identifying the third transition peak. 

For the simulations a range of out-of-subspace coupling 
strengths 7 was chosen for each of the trial Hamiltonians ( |29] l, 
(I2TI 1 and (|30] |, and the corresponding subspace leakage rate e 
as well as the analytical lower bound ej computed. For each 
of the Hamiltonians we then simulated experimental Rabi data 
and computed the discrete Fourier spectrum. The observation 
time in all cases was 30 Rabi cycles and the number of en- 
semble measurements was = 1024. The number of en- 
semble measurements for the Rabi data simulations was grad- 
ually increased until a statistically significant third peak was 
found ( |32] |. or dSTl i was satisfied, respectively. 

Fig. I2] shows the number of ensemble measurements 
necessary to conclude that the system is imperfect in the sense 
that leakage is statistically significant for the three-level sys- 
tem governed by ( [29] l for both methods. The horizontal axis 



represents the analytical value of confinement £(7). Both 
curves scale roughly l/>/lVe, which is consistent with the 
scaling of the projection noise, and hence the errors associ- 
ated with estimating e; and detecting a statistically significant 
third peak. For the three-level system it is clear that confirm- 
ing imperfect confinement by verifying dSTI ) requires more en- 
semble measurements than detecting a third peak according 
to ( |32] ). This is not too surprising since for a three-level sys- 
tem there is only one additional transition |0) ^ |2), and from 
the derivations of the confinement equations (|9]l we have. 



a=0 
N-1 



(33) 



a=0 



i.e., there is a conservation law for the cumulative sum of all 
the peak heights. Hence, if the number of possible additional 
peaks is small, then for a given level of confinement, the ad- 
ditional peaks will be greater, and thus easier to detect, than 
for a system with weak coupling to a large number of out- 
of-subspace levels, and hence many small transition peaks. 
We therefore conjecture that estimating subspace leakage us- 
ing dSTT i will become preferable for a system with coupling to 
multiple out-of-subspace levels. The results of numerical sim- 
ulations for the Hamiltonians (I2TI ) and ( l30l l. shown in Fig [8] 
support this conjecture. We observe the same general scal- 
ing behavior as for the three-level system. For the four-level 
system it is clear that although searching for the additional 
transition peak is still somewhat more efficient, the difference 
between both methods is small. For the six-level the curves 
have swapped position, i.e., using the confinement equations 
has become a more efficient way to ascertain statistically sig- 
nificant subspace leakage. 

In Appendix A we have included simulations for similar 
Hamiltonians up to ten levels to show the effective crossover 
of the curves and how the efficiency difference between the 
two methods increases with the number of additional levels. 
Note that for all the simulations we have endeavored to look 
at approximately the same range of subspace leakage. From 
these simulations it is clear that searching for the third peak 
in the Fourier spectrum is only really beneficial for systems 
with at most one extra transition. Hence, the proposed method 
for estimating subspace leakage will be more efficient than 
obvious alternatives in most cases. 



VI. THE EFFECT OF DECOHERENCE 

It is well known that even if subspace leakage is theoret- 
ically suppressed for an arbitrary control field, it is unlikely 
that decoherence will also be suppressed. Hence, we need 
to examine if the proposed confinement protocol will still be 
effective in the open system case when a qubit is subject to 
decoherence, possibly of the same order, or greater, than sub- 
space leakage. 
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a) b) 




0.0025 0.0055 0.0097 0.0113 0.0147 0.002 0.006 0.0076 0.011 0.014 

Confinement (f ) Confinement (e) 

FIG. 8: Number of ensemble measurements required to ascertain significant subspace leakage (imperfect confinement) for the four-level 
system governed by | |2U [Fig. a] and the six-level system governed by l |30| l [Fig. b] using the confinement equations and identifying a third 
peak. 



The study of arbitrary decoherence for A^-level systems is 
a lengthy discussion, including Markovian and possible non- 
Markovian processes. Even for the simpler case of Markovian 
decoherence we would need to consider the complete iV-level 
decoherence model with all the associated restrictions of com- 
pletely positive maps ll32ll . Hence, we will instead only focus 
on a restricted case to show that, for a simple example, deco- 
herence does not invalidate the protocol. It should be stressed 
that this only represents a preliminary analysis under a spe- 
cific model of decoherence. Further work will involve investi- 
gating more complicated and system-specific decoherence ef- 
fects such as A^-level dephasing and spontaneous emission as 
well as possible system specific non-Markovian decoherence. 
However, due to the extremely complicated nature of such an 
analysis we will limit our discussion to a specific case. 

We consider a perfectly confined qubit which undergoes 
Markovian decoherence and hence can be described by the 
quantum Liouville equation 

3 

dtp^~UH,p]+Y,'CkCM (34) 

where, Ck[p\ = {[Lk-, pL\] + [LkP, L\])/2, H represents the 
single qubit control Hamiltonian, and Lk are the Lindblad 
quantum jump operators, which describe the effect of the en- 
vironment on the system, each parameterized by some rate 

rfc>o. 

For a basic decoherence analysis we restrict the Lindblad 
operators to the Pauli set, {Lk] — {X^ F, Z}, and consider a 



perfectly confined, control Hamiltonian of the form 



H = '^[cos{9)Z + sin{e)X]. (35) 

This decoherence model is sufficient to describe pure dephas- 
ing as well as symmetric population relaxation processes in 
any basis, although not asymmetric relaxation processes. In- 
cluding each Pauli Lindblad term with an associated decoher- 
ence rate eliminates the problem of a preferential basis for 
qubit decoherence since any basis change of the overall sys- 
tem will only act to change the form of the Hamiltonian. 

We can solve the master equation under this model by using 
the Bloch vector formalism. Expressing the density matrix as 
p{t) = I /2 + x{t)X + y{t)Y + z{t)Z, Eq. ^ takes the form 
dtS{t) = AS{t), where S{t) = {x{t), y{t),z{t))'^ and 



/-2(r,y + r,) -dcos(0) \ 
A=\ dcos{e) -2(r^ + r,-) -dsin{e) \. 

\ dsin(6') -2{T^ + Ty)) 

(36) 

The Rabi oscillations under this evolution are described by 
the function f{t) = Tr[Po/3(t)] = (1/2) + z{t), where 
Po = |0)(0|, with an initial state p(0) = |0)(0| =i> 
5(0) = (0, 0, 1/2)^. Taking the Fourier transform of f{t) 
leads to the rather complicated general expression ( IB7l i in Ap- 
pendix B. The real component of this function describes three 
Lorentzians centered about a; = and uo = ±d. Assuming 
that d ^ Fj; J, 2, we can expand Eq. ( IB7I ) around a; = and 
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±d to obtain the functions [See Appendix B], 



/^o,l = 



81112(6*) 



cos2(6i) 



2 w2 

r/3 



rr 



(37) 



4 ± rf)2 + r2 ■■ 



where r„ = 2(ry + T,+ coa^{e){T, - T,)), = r,(l + 
sin2(6')) + ry +r3(2-sin2(6')) and ho contains a (5(w) offset 
due to the fact we are measuring the observable Pq. In order 
to describe how the maximum peak of each Lorentzian varies 
with r we integrate ho and ho^i around an interval ?y of the 
peak height 

^ , , cos2(0) p ^ 1 r ^ N 



(38) 



1 cos2(e') / ?7 
= - H arctan — , , 

2 TT \i-a/ 



sin- (61) / T] 



Hence, under decoherence the peak heights in the Fourier 
spectrum vary as a function of the integration window 77 
and the decoherence rates Ta^p- This is consistent since 
as Ta^p 0, both arctan functions approach 7r/2 and 
Hq + 2hQ i — 1. The integration window 7] is analogous to 
frequency resolution of the Fourier transform Auj, while the 
total area of the Lorentzian is equal to the peak heights when 

r — n 



Hence for small T^^y.z we can simply choose 
the resolution of the Fourier transform such that the entire 
Lorentzian is essentially contained within the data channel of 
the primary peak. 

Consider the case where we wish to ensure that the sub- 
space leakage does not exceed (. Using the upper bound for 
the subspace leakage dTSI l we have, assuming that the integra- 
tion interval is approximately equal to the frequency resolu- 
tion of the DFT (i.e. 7] w Auj) 



^2ha{Auj) + AhQ,i{Auj) 



1 



(1 - 2C)2 + 1 



41-20^ 



cos2(6i) fAuj 
arctan 

TT \Ta 

arctan 



arctan 



Auj 



2 V r 

(39) 

Here the last line assumes that Ta ~ = T. When the Rabi 
frequency is much greater than the inverse of the decoherence 
rate (as necessary for any qubit realistically considered for 
quantum information processing), then the entire Lorentzian 
broadening caused by decoherence will be contained within 
one frequency channel. Thus, Eq. (|39] l allows us to calculate 



the maximum frequency resolution of the Fourier transform 
for successful leakage estimation using our protocol. For ex- 
ample, if r w lO^^'^s^^ and we wish to confirm that the sub- 
space leakage is at most e^ax = 10^^, then the resolution of 
the Fourier transform cannot exceed A/ « 250Hz if only the 
primary peak channels are used. Obviously, this restriction on 
the frequency resolution can be lifted by including multiple 
channels around the central peak when estimating the peak 
area. 

Although the decoherence model considered is not the 
most general possible case for an imperfectly confined con- 
trol Hamiltonian, this calculation demonstrates that the ef- 
fect of decoherence does not void the protocol for estimat- 
ing subspace leakage for a common decoherence model. A 
more detailed analysis considering a full iV-level decoherence 
model, including the effect of spontaneous emission and ab- 
sorption processes and the possibility of system-specific non- 
Markovian decoherence is desirable but beyond the scope of 
the current paper. 



VII. CONCLUSIONS 

We have introduced an intrinsic protocol for "quantify- 
ing" the degree of subspace leakage for a realistic 'qubit' 
system. The protocol relies on very minimal theoretical as- 
sumptions regarding qubit structure and control, and utilizes 
a measurement model that is restrictive but extremely com- 
mon to a wide range of qubit systems. We have introduced 
a quantitative measure of subspace leakage, and shown that 
the discretization noise as a result of finite sampling does not 
limit the ability of the protocol to quantify (with appropri- 
ate error/confidence bounds) the subspace leakage for well- 
confined (near perfect) qubits. 

The ability to experimentally characterize subspace leak- 
age to a high degree of accuracy using automated, system in- 
dependent methods, which rely on the intrinsic control and 
measurement apparatus of the quantum device (required for 
standard quantum information processing) will be vital for the 
commercial success of quantum nano-technology. This pro- 
tocol represents one of the first steps in a general library of 
characterization techniques that will be required as "quality 
control" protocols once mass manufacturing of qubit systems 
becomes common. 

Although, in this discussion, the qubit state |1) is only de- 
fined through the strongest transition it should be emphasized 
that if confinement estimates are made on multiple control 
fields (for example two separate Hamiltonians which induce 
orthogonal axis rotations), the computational |1) state must 
be common for both Hamiltonians. This is not a significant 
problem, since for well engineered qubits, the computational 
1 1) state will be known on theoretical grounds. 

There are many open problems including subspace leakage 
estimates for systems undergoing a whole range of potential 
decoherence processes, quantifying confinement for multi- 
qubit control Hamiltonians and combining these schemes with 
other proposed methods for system characterization. Hope- 
fully, in the near future, a complete set of characterization 
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protocols will be developed which will augment large scale 
manufacturing techniques, allowing for efficient and speedy 
transition of quantum technology from the physics laboratory 
to the commercial sector. 



APPENDIX A: EFFICIENCY COMPARISON FOR 
LEAKAGE DETECTION PROTOCOLS 

The following simulations examined the minimal number 
of ensemble measurements required to detect imperfect qubits 
either via the confinement equations or by directly detecting 
the third transition peak. Three-level, four-level and six-level 
Hamiltonians are found in the main text, the additional simu- 
lations were performed for all other multi-level systems up to 



ten levels. The general form of each of the trial Hamiltonians 
are subsets of the ten-level system, 

9 

Hw = Y.Ek\k){k\+-ikm{k\ + |fc)(0|) (Al) 

where {Ek} = {0, 1, 1.5, 1.7, 1.9, 2.2, 2.5, 2.7, 3, 3.2}, 71 = 
1 and 7fc = 7 for fc 7^ 1. 

For each lower level system the appropriate Hamiltonian is 
simply formed by removing the appropriate number of rows 
and columns from Hiq (i.e. compare H^^ and Hq in Eqs. ( 1211 1 
and (I30I1). Each of these systems were simulated leading to 
the following results [Figs l9l [TOl and [TTl . 




FIG. 9: Number of ensemble measurements required to ascertain significant subspace leakage (imperfect confinement) for the five-level system 
[Fig. a] and the eight-level system [Fig. b] using the confinement equations and identifying a third peak. 




FIG. 10: Number of ensemble measurements required to ascertain significant subspace leakage (imperfect confinement) for the seven-level 
system [Fig. a] and the nine-level system [Fig. b] using the confinement equations and identifying a third peak. 
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U sing confinement equations 
Identifying a tliird peak. 




.0051 .0078 

Confinement (e) 



.0111 .0148 



FIG. 1 1 : Number of ensemble measurements required to ascertain 
significant subspace leakage (imperfect confinement) for the ten- 
level system using the confinement equations and identifying a third 
peak. 



noting that 



C[f{t);Lo]~2S[f{t);io]= / f{t)e'^-'^T+[fit);Lo]. 

Jo 

(B2) 

Taking the sine and cosine transforms of dtS{t) = AS{t), 
noting that 



C[f{t);uj]^ujS[f{t);uj]-f{0), 



ujS[S{t);Lu] - S{0) = AC[S{t);w], 
~ujC[S{t);uj]^ AS[S{t)-w]. 



(B3) 



(B4) 



APPENDIX B: SOLUTIONS TO THE DECOHERENCE 
MASTER EQUATION 

Here we show the derivations of Eq. [37] by solving the 
Bloch equation dtS{t) = AS{t), with A given in Eq. ( |36l ). 
To solve this differential equation, we convert to Fourier 
space. Since the Fourier transform for a system governed by 
decoherence-induced semi-group dynamics is only defined for 
t > 0, we use the cosine and sine transforms 



= / f{t)cos{iot), 

Jo 

J poo 
/(i)sinM), 
9 



(Bl) 



Combining these equations and setting S{uj) — T+[S{t)\uj] 
yields. 



and hence 



iujS{uj) - 5(0) = AS{oj), 



S{lu) = -{A-iujI)-^S{0). 
The initial condition S{0) = (0, 0, 1/2)^ thus gives. 



(B5) 



(B6) 



FT[zit)] 



c'^<f + (2r^ + 2r^ + iuj){2Ty + 2r^ + iuj) 



2{c'^<P + (2r^ + 2r^ + iuj){2Ty + 2r^ + iu;))(-2(r^ + Vy) - ioj) - 2(Ps^{2Ty + 2T^ + iuj) ' 

I 



(B7) 



where c = cos{d) and s = sin(6'). The subsequent expan- 
sions are too lengthy to include here, however standard sym- 
bolic toolkits such as Mathematica can handle such expres- 
sions. The first step is to consider only the real component of 
FT[z{t)]. Next, the denominator is expanded to second order 
around = or = ±d. After this, we expand the numera- 
tor and denominator, neglecting all terms of the form Tx,y,z/d 
and smaller, assuming T^^y^z <C d and being careful to note 
that for expansions around u; = ±d we must keep terms of the 



form ujTx^y^z/d. After simplifying the expressions we find 



COS^{0) 



ho,i = 



. F2 ' 

a 

r/3 



{uo ± df + VI ' 



(B8) 



where = 2{Ty+T ^ + cos^{e){T .^-T ^)) anAT p = r^(l + 
sm^{e)) + Vy+ r^(2 - sin^{e)). Confirming that Eq. ^ 
describes three Lorentzian curves centered on = and lo = 
±d. 



13 



ACKNOWLEDGMENTS 

SJD acknowledges the support of the Rae & Edith Bennett 
Travelling Scholarship. SGS acknowledges support from an 
EPSRC Advanced Research Fellowship and the Cambridge- 
MIT Institute. DKLO acknowledges support from Sidney 
Sussex College, Cambridge and SUPA. SGS and DKLO also 



acknowledge support from the EPSRC QIP IRC (UK). SJD, 
JHC and LCLH are supported in part by the Australian Re- 
search Council, the Australian Government and the US Na- 
tional Security Agency (NSA), Advanced Research and De- 
velopment Activity (ARDA), and the Army Research Office 
(ARO) under contract number W911NF-04- 1-0290. 



[1] J. E. Mooij, T.P. Orlando, L. Levitov, L. Tian, C.H. van der Wal 
and S. Lloyd, Science. 285, 1036-1039 (1999). 

[2] Y. Nakamura, C. D. Chen and J. S. Tsai, Nature (London). 398, 
786 (1999). 

[3] M. R. Geller, E. J. Pritchett, A. T. Sornborger, R K. Wilhelm, 

quant-ph/0603224 (2006). 
[4] M.D. Lukin and P.R. Hemmer, Phys. Rev. Lett. 84, 2818 (2000). 
[5] A.D. Greentree et. ai, J. Phys: Condens. Matter. 18, (2006), 

S825-S842. 

[6] J. Wrachtrup, S. Y. Kilin and A. P. Nizovtsev, Optics and Spec- 
troscopy. 91, 429-437 (2001). 
[7] J.I. Cirac and P Zoller, Phys. Rev. Lett. 74, 4091 (1995). 
[8] J. Clarke et ai. Science. 239, 992 (1988) 
[9] A. Blais, A. Maassen van den Brink and A. M. Zagoskin, Phys. 
Rev. Lett. 90, 127901 (2003). 
[10] A.M. Steane, Phys. Rev. Lett. 77, 793, (1996). 
[11] R. Laflamme, C. Miquel, J. P Paz, and W. H. Zurek, Phys. Rev. 

Lett. 77, 198 (1996) 

[12] D. Gottesman, Ph.D Thesis (Caltech), [quant-ph/9705052] 
[13] L.-A. Wu, M.S. Byrd and D.A. Lidar, Phys. Rev. Lett. 89, 
127901 (2002). 

[14] M.S. Byrd, D. A. Lidar, L-A. Wu and P Zanardi, Phys. Rev A. 

71, 052301 (2005). 
[15] R. Fazio, G. Palma and J. Siewert, Phys. Rev. Lett. 83, 5385 

(1999). 

[16] C. Mochor, Phys. Rev. A. 69, 032306 (2004). 
[17] J. Preskill, Introduction to Quantum Computation, pages 213- 
269. World Scientific, Singapore, 1998. 



[18] M. Grassl, Th. Beth, T. Pellizzari, Phys. Rev. A. 56, 33 (1997). 
[19] J. Vala, K.B. Whaley and D.S. Weiss, Phys. Rev. A. 72, 052318 
(2005). 

[20] K. Khodjasteh and D.A. Lidar, Phys. Rev. A. 68, 022322 

(2003) . 

[21] P Aliferis, B.M. Terhal, Quant. Inf. Comp. 7, 139-156 (2006). 

[22] J.H. Cole etal., Phys. Rev. A. 71, 062312 (2005). 

[23] J.H. Cole, S.J. Devitt and L. C.L. Hollenberg, J. Phys. A: Math. 

Gen. 39, (2006), 14649-14658. 
[24] S.J. Devitt, J.H. Cole, L.C.L. Hollenberg, Phys. Rev. A. 73, 

052317 (2006). 
[25] J.H. Cole et. al., Phys. Rev. A. 73, 062333 (2006). 
[26] M. H. Devoret and R. J. Schoelkopf, Nature (London). 406, 

1039 (2000). 

[27] V. I. Conrad, A. D. Greentree, D. N. Jamieson and L. C.L. Hol- 
lenberg, J. Comput. Theor. Nanosci. 2, 214 (2005). 

[28] A. Aassime, G. Johansson, G. Wendin, R. J. Schoelkopf, and P. 
Delsing, Phys. Rev. LeU. 86, 3376 (2001). 

[29] J. M. Martinis, S. Nam, J. Aumentado and C. Urbina, Phys. 
Rev. Lett. 89, 1 17901 (2002). 

[30] G. Wendin and V. S. Shumeiko, lcond-mat/0508729 (2005). 

[31] R. N. Bracewell, The Fourier transform and its applica- tions, 
McGraw-Hill series in electrical and computer en- gineering. 
Circuits and systems. (McGraw Hill, Boston, 2000), 3rd ed.l4. 
(Clarendon Press, Oxford, 1990). 

[32] S.G. Schirmer and A.I. Solomon, Phys. Rev. A. 70, 022107 

(2004) . 




0.0025 0.0055 0.0097 0.0113 0.0147 

Confinement (e) 



